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LEAD ISOTOPES AS A KEY TO THE RADIOACTIVITY 
OF THE EARTH’S MANTLE 


BY 
R. D. RUSSELL 


Introduction 


At present, the earth can be divided into three major parts: the core, 
the mantle, and the crust. Of these, only the crust is amenable to close 
investigation. There now seems to be good evidence, however, that the 
crustal rocks have grown during geological time by additions of volcanic 
material derived from the upper part of the mantle. If this is true, the 
crust cannot be considered to be a true primary subdivision of the earth, 
and a detailed knowledge of the properties of the mantle becomes even 
more important for a proper understanding of the very early history of the 
earth. 

Astronomical measurements of the orbits of the earth and moon give 
the mass and moment of inertia of the earth and, therefore, place limita 
tions on the distribution of mass within the mantle. The study of seismic 
waves from earthquakes have given the elastic constants and temperature 
distribution within the mantle (Jeffreys, 1952; Jacobs, 1953). Such 
studies indicate that the mantle is made preponderantly of ultrabasic 
rocks, but no way has been found of discovering the details of its chemi- 
cal composition. 

A particularly important question in this respect is the proportions 
of radioactive isotopes in the mantle. On this rests our proper under- 
standing of the entire thermal history of the earth. Without this data, 
it is impossible to decide whether the earth is cooling or heating, and 
this must be established before the various theories of mountain building 
can be properly assessed. 

Studies of the thermal history of the earth are dependent on assumed 
values for the distribution and concentration of uranium, thorium, and 
potassium in the mantle. In the past, these values have been obtained by 
analogy with certain rock types and with various classes of meteorites 
assumed to represent different parts of a destroyed planet. Urey (1952) 
has given evidence that casts doubt on this interpretation of meteorites. 

While the knowledge of radioactivity of the mantle is extremely 
important, and while very little is known about it, any means of casting 
light on the subject seems worth while. Some tentative information on 
this subject has been obtained by the writer from studies of common 
lead abundances. The results require verification that can be provided 
only by many more isotopic analyses than are available now. They are 
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presented here to illustrate the part they may be expected to play in 
studying the radioactivity of the earth’s mantle. 


The Systematics of Common Lead Isotope Abundances 


The source of lead ores. Early studies of the atomic weight of 
lead found in uranium and thorium minerals established that it differed 
greatly from that of lead found in lead minerals. Subsequent spectro- 
scopic (Rose and Stranathan, 1936) and mass spectroscopic studies 
(Aston, 1936; Mattauch and Hauk, 1937) showed the great differences 
in various radiogenic leads, and also the differences between such leads 
and common leads. It was much more recently, however, that variations 
in the isotopic abundances of common leads were discovered and studied 
in detail (Nier, 1938; Nier, Thompson, and Murphey, 1941; Collins et al., 
1953; Russell et al., 1954; Allan et al., 1953; Geiss, 1954). 

The interpretation of common lead abundances has been attempted 
at various times (e.g. Holmes, 1946, 1947; Houtermans, 1947; Alpher 
and Herman, 1951). All of these interpretations are based on the funda- 
mental assumption that all terrestrial lead, at some early stage, was 
well mixed in a hot and active earth, and had everywhere the same 
isotope ratios. They suppose further that to the primeval lead in the 
outer part of the earth was added radiogenic lead from the uranium and 
thorium with which it was associated, but that no other process operated 
within the earth to change the lead isotope ratios. When such lead was 
formed into lead minerals, it became isolated from uranium and thorium, 
and its abundances were preserved unchanged unless the minerals were 
subsequently dissolved and taken again into contact with the radioactive 
elements. Such an interpretation works so well that there can be little 
doubt that it is very close to the truth. 

Mathematical analyses of lead isotope abundances are based on 
either of two additional assumptions. One, due to Alpher and Herman 
(1951), supposes that the source rocks from which the lead ores were 
derived are approximately homogeneous in uranium, thorium, and lead. 


The other, due to Holmes (1946), assumes that uranium to lead ratios 


differ from place to place but that, in each place, the ratio has not been 
disturbed by additions of, or depletions of, either uranium or lead except 
by simple radioactive decay. Both treatments have been developed fully 
and, although neither assumption is exact, both seem to give a remarkably 
good representation of measured lead abundances. 

The fact that these two different approaches do work provides two 
important clues about the origin of lead ores. These clues are: (1) lead 
ores are derived chiefly from a source in which uranium, thorium, and 
lead are approximately uniformly distributed; and (2) the source material 


a 
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must not have been mixed on a large scale at any time since the final 
consolidation of the outer part of the earth in such a manner as to change 
local uranium to lead ratios. 

It is known that the rocks forming the continents are of a very 
heterogeneous nature, and that their radioactive contents vary greatly. 
Moreover, all radioactive age determinations strongly suggest that the 
continents have been built progressively during geological time, and that 
they are in fact still growing. Thus continental rocks do not satisfy the 
above conditions for the source rocks of lead ores. It is possible, however, 
that the upper part of the mantle is more homogeneous than the crust, 
and that it has been less disturbed during the earth’s history. Thus our 
present knowledge of the behavior of lead isotope abundances and the 
nature of the earth makes the derivation of lead ores from the upper part 
of the mantle seem plausible, and their derivation from crustal rocks 
highly improbable. Since no third alternative seems to exist, the con- 
clusion that the lead in most lead ores has been derived from the upper 
part of the earth’s mantle seems inescapable. Consequently, the study of 
lead isotope abundances should provide information about the uranium 
(and thorium) distribution within the outer part of the mantle. 

A difficulty in the simple theory. Following the method of Alpher 
and Herman (1951), but adopting the nomenclature of Collins, Russell, 
and Farquhar (1953), lead existing in the mantle at time t, and hence in 
an ore separated from the mantle at time t, should be represented in 
isotopic constitution by the formulae 


=a- 137.8V (e* - 1), 


x 
a ee V (e** — 1), 
a eats Ww (e** ~ 1), (1) 
TABLE 1 


DEFINITION OF SYMBOLS 


Pb? /Pb*°4 
Pb?°7 /Pb?°4 


Pb?°8 /Pb?°4 
y238 /Pb*4 
y735 /Pb?°4 
Th?232 /Pb?°4 


X =0.154 x10° years* 
N’=0.972 x 10° years* 
NV’ = 0.0499 x 10° years* 
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where the symbols used are defined in TABLE 1. In general, the unknown 
parameters V, W, a, 6, and c, are chosen by fitting these equations to 
the observed abundances of leads of known age, although one other 
method has been attempted (Allan et al., 1953). In all of these cases, 
the form of EQUATIONS 1 was assumed, including the assumption that 
the V appearing in the first two equations were identical. It occurred to 
the writer that a good test of the assumptions could be made by determin- 
ing separately values for V from each of the first two equations using 
abundances of the same dated lead samples, and comparing the results. 

For this purpose, samples from three main geological divisions of 
the Canadian shield were used. Several samples were chosen from each 
region in the hope that their mean age might be equated more accurately 
to the age of the corresponding orogenies. 

The first nine samples occurred in the Grenville geological province. 
The Grenville province and its relation to some other major geological 
divisions of North America are described as follows by King (1951, p. 4): 
“Younger orogenetic belts thus succeed older belts on passing outward 
from the core of the continent. This pattern is especially well dis- 
played toward the southeast, where one of the later Pre-Cambrian oro- 
genic belts, the Grenville, borders the ancient deformed rocks of the 
core and is followed farther out by the Appalachian belt, of Paleozoic 
age, which parallels the Grenville belt for part of its course.’’ 


This supposition that the Grenville represents a single orogenic 
belt is now rather widely accepted. Wilson (1954) suggests that the 
igneous activity accompanying the building of a single orogenic belt is 
confined to a period of about 200 million years duration. This is sup- 
ported by many age determinations on uraninites, feldspars, and micas 
from Grenville pegmatites and on one Grenville rock, essentially all of 
which give ages between 850 million years and 1050 million years (Ells- 
worth, 1932; Nier, 1938; Nier, Thompson, and Murphey, 1941; Collins 
et al., 1954; Wasserburg and Hayden, 1955; Shillibeer and Watson, 
1955; Tilton et al., 1955; Cumming et al., 1955). Leads found in 
such an area could conceivably be much older or much younger than the 
age of pegmatites and of the country rocks, but studies of common lead 
abundances have been consistent with the supposition that the time of 
formation of most lead minerals coincides with the earlier stages of 
orogenic activity. Thus an age of 1000 million years has been assigned 
to the nine Grenville leads. 

“The region north of Lake Superior—mainly in the province of 
Ontario, but extending northeastward into Quebec and northwestward 
into Minnesota and Manitoba—is formed of a complex of disturbed and 
metamorphosed sediments and volcanics, and of granitic plutonic rocks. 
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These are probably the most ancient rocks in the Laurentian Shield’? 
(King, 1951, p. 6). Eight isotopically analyzed leads were chosen from 
various parts of this area, which is called by Wilson (1954) the 
Keewatin province or continental nucleus. They were assumed to have a 
mean age of 2500 million years. This age was obtained by averaging two 
lead ratio ages of 2470 and 2540 million years found by Nier (1939; and 
Nier, Thompson, and Murphey, 1941) for a monazite and a uraninite from 
the Huron Claim on the Winnipeg River in southeast Manitoba. A second 
uraninite, collected by the writer from the same pegmatite gave a lead 
tatio age of 2550 million years in excellent agreement with Nier’s re- 
sults. This age is also supported by potassium-argon ages found by 
Shillibeer and Russell (1954). 

Although it is admittedly a big assumption to assume that this age 
applies to the whole Keewatin, scarcity of uranium minerals throughout 
this geological province have prevented more detailed studies of its age. 
Since the oldest lead-ratio ages reported for minerals anywhere are 2650 
million years (Holmes, 1954), it seems unlikely that the age of 2500 
million years assigned to the Keewatin has been underestimated greatly. 

The Yellowknife province in Canada’s Northwest Territories is 
similar petrologically to the Keewatin, and it is also believed to be a 
continental nucleus (Wilson, 1954). No lead ratio ages have been re- 
ported for this area because of the extreme rarity of Yellowknife uranium 
minerals. A single potassium-argon measurement on the Moose Dyke 
pegmatite, however, indicates an age of 2150 million years. While dif- 
ferent methods for dating common leads sometimes yield very different 
results, all agree in giving ages for Yellowknife galenas less than those 
for Keewatin galenas (Bullard and Wilson, 1954). Thus an age of 2100 
million years for the Yellowknife leads seems to be reasonable. The 
effect of adjusting the ages assigned to these three groups of leads will 
be discussed below. 

The 28 lead samples for these regions, with their Pb2°°/Pb 24 and 
Pb 207/Pb 2% ratios are listed in TABLE 2. The isotope ratios, with the 
assigned ages, were used to determine the parameters in the equations: 


x =a- 137.8 V (e* -1) 

eo Voie 4), (2) 
The constants V and V’ can be found by plotting x against 137.8 (ert - 1) 
and y against (er* — 1), and finding the negative of the slopes of the 


straight lines that should result. This was done by the methed of least 
squares with the results: 


V =0.075 +0,005 
V’=0.096 + 0.005 
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In calculating V and V’, it was assumed that all the uncertainty in the 
experimental numbers was associated with the ages t of the ores, and 
that the abundances x and y were exact. This should be a good approxi- 
mation. The limits of error indicated are standard deviations. 

The two values V and V’, which previously always have been taken 
as identical, differ by more than four times the standard deviation of 
either. This difference appears to be real. Four possible explanations 

for this discrepancy are discussed below: 


(1) One or more of the constants used may be in error. Three are 
used in the calculations, two of which would have to be in error by an 
impossibly large amount to explain the effect. The first of these, the 
ratio of the uranium isotopes (137.8), has been determined several times 
and no values differing from the one quoted by more than 1 per cent have 
been reported. An error of 25 per cent would be required to explain the 
discrepancy. An error in the decay constant for U*® of 20 per cent 
could also explain the discrepancy, but many recent measurements 
indicate a precision of a small fraction of 1 per cent. An error of only 
8 per cent in the decay constant for U*> would explain the difference, 
while the accepted uncertainty of the value used is between 1 and 2 per 
cent. Raising the decay constant by this amount, however, would de- 
crease the ages assigned, which were based primarily on the lead-ratio 
method. As shown below, younger ages would aggravate the difficulty. 

The difference cannot be explained on the basis of mass spectrometer 
errors which would occur primarily in the measurement of the Pb 2% 
abundances and, consequently, change V and V’ proportionately. 


(2) The ages assigned to the samples might be in error. The nature 
of the equations are such that changing the age assigned to the Grenville 
samples changes V more rapidly than V’, while changing the age for the 
Yellowknife and Keewatin samples changes V’ more rapidly than V. Thus 
the difference between the two calculated values is diminished by in- 
creasing the ages assigned to any or all of the samples, and is increased 
by decreasing the ages assigned to any or all samples. 

If the Grenville age alone is in error, it must be low by about 30 per 
cent. Since the Grenville is undoubtedly the best dated Pre-Cambrian 
province anywhere, its age being determined by dating many pegmatites 
and one granite, such an error seems unlikely. No isotopic lead or 
potassium-argon ages greater than 1100 million years have ever been 
reported for this region. Although it is possible that Grenville leads may 
“be older than the Grenville orogeny, it seems unlikely that nine samples 
selected at random should all have been preserved from an earlier 
orogeny. It seems more likely that the age assigned, if in error, is high. 

Greater uncertainty is attached to the Keewatin-age. This age is 
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based on the lead-ratio method which the writer believes preferable to 
all others for dating Pre-Cambrian minerals. However, the lead-ratio 
method has been attacked strongly at the same time by others, who con- 
tend that sufficient radon may be lost from uranium minerals during their 
history to cause an appreciable depletion of Pb 7 (Kulp et al., 1954). 
If this latter argument were to be accepted, the assumed age for the 
Keewatin should be lowered, and the disagreement between V and V’ 
increased. To the writer’s knowledge, no one has suggested that the 
lead-ratio method gives too low a result. 

Since the Keewatin is believed to be a continental nucleus, the 
earth being presumed too unstable to preserve continents prior to that 
time, there seems little likelihood that the Keewatin galenas are signif- 
icantly older than Keewatin rocks. 

Arguments similar to those for the Keewatin apply to the Yellow- 
knife. The age, which is calculated by the potassium-argon method with 
a branching ratio of K* of 0.089, would be lowered if the branching ratio 
were to be increased. Although a higher branching ratio is preferred by 
many, there is no indication that a lower value would be better. Again, 
the Yellowknife province is a continental nucleus and, consequently, 
should not contain minerals from some earlier epoch. 

Thus it seems unlikely that the ages of these three geological 
provinces have been underestimated, and it is even more unlikely that 
the mean age of the randomly chosen galenas should be appreciably 
older than the age of the surrounding country rock. 

(3) The mantle, like the crust, may have been changing during 
geological time, and chemical and physical processes other than radio- 
activity may have changed its uranium content relative to its lead con- 


tent. This would invalidate Equations 1, which should then be re- 
placed by the equations: re 
x=a-137.8/f Vit)e dAdt, 


pid seth id. wenftcVit era Coe (3) 


EQuaTions 3 become identical to the simple EquaTions 1 if V(t) is 
assumed constant. Remembering that the values of V and V’ obtained 
before were obtained from the slopes of the x versus 137.8 (ert — 1) and 
y versus (e* * — 1) plots, we find that for the more general EQquaTIONS 3 


ws —dk Bhd ek gente. 


Sete task Se Vit) 
137.8 d(e“1) — 137.8re**_ dt 


pe ee er wa dy 
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Applying these formulae to the more general EquaTions 3 then shows 
that even if V(t) is not constant, its value should be given by the slope 
of either of the two plots mentioned above. The slopes calculated by 
least squares would then both represent averages of V(t) over the range 
of time considered. Thus these averages should be essentially the same, 
and the difference between V and V’ cannot be explained by assuming any 
arbitrary depletions or enrichments of uranium or lead in the mantle 
during the interval of time between the formation of Keewatin and Gren- 
ville galenas. 

(4) A final possibility is that the older and younger leads were de- 
rived from different parts of the mantle that had essentially different 
uranium to lead ratios preserved at each place throughout the history 
of the earth. 

The effect of this assumption can be seen approximately without 
working through the details of the calculations by noting that the old 
leads (grouping the Keewatin leads and Yellowknife leads together for 
convenience) would then satisfy the equations 


Xx. = a, + 137.8 V, (eAto — erty, 


1 ° 


and ¥, = b, + LA (e' ~et1), 


oO 
and the Grenville leads 


x, = a + 137.8V, (eMo —e4*2), 


Bits ely pe bo. + V, (eAto ~e*%), (5) 
Plotting the abundances of leads satisfying EQUATIONS 5 in the manner 


previously described, and calculating the negatives of the slopes as 
before, one would obtain values related to V, and V, by the equations 


pe ee tae eer wa a ae 
137.8(e1 — et2) 


(V, — V; Yerto tr A M1 = V, eM2 
(M41 - e M2) ¥ 


? 


(V, =i, Jerto+ V, eM 4 =<K, ehty (6) 


and Y= - ; 

(At, — At) 
Unless V, is equal to %, V will differ from V’ Substituting the 
values obtained above for % and V’, and using for t, the value of 4200 
million years given by Russell and Allan (in press), values for and V2 
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were determined from EQuaATiIoNs 6 and found to be 


V; 0.062, 
and V,> = 0.067. 


These correspond to present lead to uranium ratios of 8.6 for the part of 
the mantle from which the Keewatin and Yellowknife leads were derived, 
and 8.0 for the part of the mantle from which the Grenville leads were 
derived. 

This explanation seems to be more likely than the other three. 


Conclusions 


Values tor the present ratio U5/Pb*™* for dated common minerals 
calculated from the Pb™/Pb™™ and Pb*®’/Pb™ ratios differ by far 
more than can be attributed to experimental error. This difference ap- 
pears to be real and suggests some major fault in the ordinary inter- 
pretation of lead isotope abundances. 

The four possible explanations considered above are the only 
possibilities that might explain the 25 per cent difference between the 
uranium-lead ratios, without at the same time admitting a complete failure 
of the basic equations used. The latter possibility seems farfetched in 
view of the considerable success with which these equations have been 
used to derive values for the age of the earth and to date lead ores. 

Of the various possible explanations, the simplest and most probable 
is that which assumes a derivation of the younger and the older leads 
from different parts of the mantle. It is suggested that these parts have 
maintained since an early time in the earth’s history permanent differences 
in their uranium to lead ratios. This explanation is in complete accord 
with other observations about lead ore abundances. Lead isotope studies 
have indicated that leads are derived from rocks approximately uniform 
in their uranium, thorium, and lead contents. The eight per cent dif- 
ference in the uranium to lead ratios for the parent rock from which the 
Keewatin and Grenville leads were derived is small enough to satisfy 
this requirement. In order that Holmes’ method for calculating the age 
of the earth yield sensible results, as it appears to do, the parent 


material for lead ores must be capable of preserving permanent hetero- . 


geneities in its uranium to lead ratios. This is precisely the hypothesis 
that explains most easily the abundances of these Canadian leads. 

If this interpretation of the lead abundances is true, it provides 
additional information about the earth’s mantle. First, the lead to uranium 
ratio of the mantle is approximately 8.3, which is significantly higher 
than most previous estimates (cf. Collins et al., 1953, who obtained 7.4). 
Second, the part of the mantle which gives rise to the lead ores has 


PE ——_ -_ 


Russell: Lead Isotopes and Radioactivity 447 


been able to preserve small heterogeneities in its lead and uranium 
concentrations which should have been destroyed if convection currents 
exist in that part of the mantle. Deep earthquake epicenters extend to 
depths of 700 kilometers, and it is often assumed that the material for 
the continents has been derived from depths this great. Lead ores appear 
to have been formed at the same time as the periods of mountain building. 
If they are brought up from the mantle from similar depths at the same 
time as the crustal rocks, the evidence of the common leads would seem 
to refute the possibility of convection currents in this part of the mantle. 

The writer is keenly aware of the difficulties in assigning correct 
geological ages and, in particular, of the difficulties in assigning ages 
to galenas. For this reason, he considers the results and conclusions 
reached here as tentative only. The calculations, however, do indicate 
clearly a major inconsistency in the present simple interpretation of 
lead abundances or in the ages assigned to the lead ores. Whether the 
explanation put forward here is correct, or whether it will be displaced 
by a better one, the final solution cannot help but add to our understand- 
ing of the history of lead ores and the radioactivity of the mantle. 
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